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Abstract 



For a large class of nonlinear evolution PDEs, and more generally, of 
nonlinear semigroups, as well as their approximating numerical meth- 
ods, two rather natural stability type convergence conditions are given, 
one being necessary, while the other is sufficient. The gap between 
these two stability conditions is analyzed, thus leading to a general 
nonlinear equivalence between stability and convergence. 



1. The General Setup 

The study of linear and nonlinear evolution systems of PDEs, with 
possibly associated initial and/or boundary problems can, as is well 
known, be dealt with in the more general framework of semigroups de- 
pending on a continuous parameter which represents time. Here, the 
study of the nonlinear equivalence for numerical methods approximat- 
ing exact solutions of nonlinear evolution PDEs will be dealt with in 
this more general framework, namely, of numerical methods approxi- 
mating rather general nonlinear semigroups. 

Definition 1.1. 

Given a normed vector space {X, \ \ ||). By a nonlinear semigroup on 
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X we mean any family of mappings 

(1.1) E{t):Xt — >X, te[0,oo) 
where 

(1.2) Xt+, CXtQX, t,se [0, oo) 

(1.3) E{0) = idxo 

while for every t,s e [0, oo) we have the commutative diagram 



E{t + s) 




□ 

An immediate consequence of (1.4) is that (1.1) takes the stronger form 

(1.5) E{t) -.Xt^Xo, te [0, oo) 

The nonlinear semigroup (1.1) - (1.5) is called well posed, if and only 
if the mapping 

(1.6) Z3{t,u)^ — >E{t)ueX 
is continuous, where 

(1-7) Z = Ue[o,oo) {t} X Xt 
Remcirk 1.1. 

The reason for not asking that Xt — X, for t e [0, oo), is to allow for 
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"finite blow up time" which is a frequent phenomenon in the nonlin- 
illustrated by a simple example like the semigroup E[t) 
of solutions of the ODE 

U'(t) = C/2(t), U(0) ^u, te [0, oo) 

in which case we can take (X, || ||) = R and Xt = {—oo,l/t). for 
t G (0, oo), while Xq = X, since for -u G we have U {t) = u/ (1 — ut). 
In this case E(t)u = U{t), for t G [0, oo). 

2. Three Versions of Stability 

The reason stability conditions for numerical methods solving PDEs 
are so important from practical point of view is very simple : such con- 
ditions are formulated exclusively in terms of the respective numerical 

methods, thus their verification does not in any way involve any knowl- 
edge of the exact solutions of the respective PDEs, and instead, it can 
be performed based alone on the given numerical methods. 

The classical example in this regard is the Lax-Richtmyer Linear Equiv- 
alence Theorem, [Lax-Richtmyer] , seen as the fundamental theorem of 
linear numerical analysis, and which says that the stability of a nu- 
merical method is both necessary and sufficient for its convergence to 
the exact solution of the respective linear evolution PDE, provided 
that the PDE is well posed, and the numerical method is consistent 
with that PDE. 

Here of course, the well posedness of the PDE docs not depend on 
any numerical method, being but an intrinsic property of the PDE 
under consideration. As for the consistency of the numerical method 
with the given PDE, that can usually be established based on suitable 
smoothness assumptions on the exact solutions of the PDE, without 
any more specific knowledge of such exact solutions. 
And then, once well posedness and consistency have been established, 
the crucial issue of convergence of the numerical method to the un- 
known exact solution of the given PDE is simply identical with the 
issue of the stability of that numerical method, stability which, as 
mentioned, is an intrinsic property of the numerical method, thus it 
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can be established without any further knowledge about the respective 
PDE or about its exact solutions. 

One of the main troubles with the Lax-Richtmyer Linear Equivalence 
Theorem, Rosinger [1-9], is that its proof depends essentially on the 
Principle of Uniform Boundedness of Linear Operators in Banach 
Spaces, a principle which simply cannot have a general enough non- 
linear correspondent. Therefore, general enough nonlinear extensions 
of the Lax-Richtmyer Linear Equivalence Theorem are hard to come 

by. 

Here however, it is important to note that in the proof of the Lax- 
Richtmyer Linear Equivalence Theorem, the implication "stable =^ 
convergent" is trivial, and the hard part is the proof of the converse 
implication " convergent =^ stable" , where the mentioned essentially 
linear principle is made use of. 

Therefore, one may expect that in any proper and wide enough nonlin- 
ear generalization of the Lax-Richtmyer Linear Equivalence Theorem 
the difficult part would be the implication "convergent =>■ stable". 
In other words, finding necessary conditions of nonlinear stability may 
appear to be the main challenge. 

Here, following Rosinger [2], we present three general nonlinear stabil- 
ity conditions which, as seen in the sequel, are closely related. 

Definition 2.1. 

Given a well posed nonlinear semigroup (1.1) - (1-7), the numerical 
methods associated with it will be continuous mappings 

(2. 1) [0, oo)xX3 {At, u) I — ^ CAtu e X 
which satisfy 

(2.2) Co = tdx 

□ 

In the sequel, we need the following auxiliary notion related to (1.7). 
A subset Z ' C Z is called regular, if and only if the following three 
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conditions are satisfied 

(2.3) Z' is closed 

(2.4) i?(t)X;+,CX;, t,sG[0,oo) 

(2.5) Ca*^ 

t+At ^ ^ t) At e [0, oo) 
where we denoted 

(2.6) X\^{u^Xt\{t,u)^Z'}, ie[0,oo) 

When it comes to stabihty, another of the main troubles with the 
Lax-Richtmyer Linear Equivalence Theorem, Rosinger [1-9], is that 
it is based essentially on the concepts of boundedness and complete- 
ness. Indeed, in the proof of the implication "convergent ^> stable", 
the Principle of Uniform Boundedness of Linear Operators in Banach 
Spaces is critically used, and as well known, that principle is not valid 
in arbitrary normed spaces, but only in Banach spaces, that is, normed 
spaces which are complete. 

As argued, however, in Rosinger [1-9], see also section 5 below, sta- 
bility, as well as convergence of numerical methods, is in fact related 
with continuity and compactness, rather than boundedness and com- 
pletness. And this fact becomes obvious starting with the next three 
definitions, and then, with the corresponding results related to the 
equivalence between convergence and stability. 

Definition 2.3. ( Local Stability ) 

A nonlinear numerical method (2.1) - (2.2) is called locally stable, if 
and only if 
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V regular C Z, T e (0, oo), compact K C X : 
3 L',p'>0 : 

(2.7) 

V At > 0, n e N, u,v e K n X '^^t : 
nAt<T, \\v — u\\<p' 

Definition 2.4. ( Distant Stability ) 

A nonlinear numerical method (2.1) - (2.2) is called distantly stable, if 
and only if 

V regular Z ' <Z Z, T G (0, oo), compact K <Z X, p > : 
3 L">0 : 

(2.8) 

V At > 0, n e N, u,v e K n X '^^^ : 

nAt <T, \\v -u\\> p =^ WCl^v - Cl^uW < L"\\v -u 

Remcirk 2.1. 

Two facts are worth noting here : 

1) There is a significant difference between local stability and distant 
stability as defined above. Indeed, in the former, nearby vectors u and 
V are involved, since they are subjected to the condition | — ii| | < p' 
which appears in (2.7). On the other hand, in the latter concept of 
stability, the vectors u and v involved are supposed to be sufficiently 
far from one another, as in (2.8) they are required to satisfy the con- 
dition ||v — > p. 

2) The above concept of distant stability is not the same with a concept 
of global stability. Indeed, as pointed out above, the vectors u and v 
involved in (2.8) cannot be arbitrary, since they are restricted by the 
condition | |i> — ii| | > p. 

□ 



u\ 
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Finally 

Definition 2.5. ( Stability ) 

A nonlinear numerical method (2.1) - (2.2) is called stable, if and only if 

V regular Z' C Z, T G (0, oo), compact K <Z X : 

3 L >0 : 

(2.9) 

V At > 0, n e N, u,v e K n X '^^^ : 
nAt<T =^ \\Cl^v - Cli^uW < L\\v - u\\ 

The relation between these three concept of stabihty is given by the 
easy to prove 

Proposition 2.1. ( Stability = Local + Distant Stability ) 

A nonlinear numerical method (2.1) - (2.2) stable, if and only if it is 
both locally stable and distantly stable. 

Remark 2.2. 

The need for considering regular subsets Z' C Z in the above three 
concepts of stability is due to the fact that even in simple cases, the set 
Z in (1.7) fails to be closed, as for instance in the example in Remark 
1.1., where 

Z = {{t,u) e[0,oo) xR I tu<l} 

3. Stability Characterization of the Convergence of 
Numerical Methods 

Definition 3.1. 
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Given a well posed nonlinear semigroup (1.1) - (1-7). A numerical 
methods (2.1) - (2.2) is called convergent to that nonlinear semigroup, 
if and only if 

V regular Z' CZ, T e (0, oo), compact K C X, e > : 
3 e>0 : 

(3.1) 

y te[o,T], At>0, ueN, ue Knx'^nx'^^^ : 
At,\t-nAt\<€ =^ \\E{t)u - Cl^uW < e 

Theorem 3.1. ( Convergence Distant Stability ) 

Given a well posed nonlinear semigroup (1.1) - (1.7) and a convergent 
numerical method (2.1) - (2.2) associated with it. Then that numeri- 
cal method is distantly stable. 

Proof. 

Assume that, on the contrary, there exists a regular subset Z' C 

Z, T ^ (0, cxo), a compact K C X and p > 0, such that for every 
j G N, there exists Atj > 0, uj e N, with UjAtj < T, as well as 
Uj, Vj e KnX ; .^j. which satisfy 

(3.2) \\CZv, ^ C'^,^u,\\ > j\\v - u\\ 

(3.3) \\vj-Uj\\>p 
However, we note that the set 

(3.4) H^{CltU I At e[0,T], neN, uAt <T, ue KnX'^^^} 
is precompact. Indeed, given e > 0, then (3.1) implies that the set 

(3.5) H, = { Cl^u I AtE[0,e], neN, nAt<T, ueKDX } 
is contained in an e-neighbourhood of 
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(3.6) K'^ {E{t)u I t e[o,T], ue Knx[} 
But 

(3.7) X ' is compact 

since in view of (2.6), it is the image of the compact set Z'n{[0, T] xK) 
under the continuous mapping (1.6). 

Further, for every given n e N, the set 

(3.8) H'^ = {Cl,u I Ate[0,T], ueK} 

is compact, since it is the image of the compact [0,T] x K under n 
iterates of the continuous mapping (2.1) - (2.2). 

We denote now 

(3.9) H', = [jH', 

where the union is taken over all n e N for which n9 < T. And then 
obviously 

(3.10) HCH.UH'^ 

thus the relations (3.5) - (3.10) complete the proof of (3.4). 

However, (3.4) implies that H is bounded, therefore, the relations (3.2) 
and (3.3) contradict one another. 

□ 

In a certain sense, a converse of Theorem 3.1. is given below by The- 
orem 3.2., for which we first need 

Definition 3.2. 

Given a nonlinear semigroup (1.1) - (1.7). A numerical methods (2.1) 



- (2.2) is called consistent with that nonlinear semigroup, if and only if 

V regular Z ' C Z, T G (0, oo), compact K C X, e > : 

3 5>0 : 

(3.11) 

Vte[o,r], At>o, ue Knx[^^^ : 

At<S =^ \\CAtE{t)u - E{At)E{t)u\ \ < eAt 

Theorem 3.2. ( Local Stability =^ Convergence ) 

Given a well posed nonlinear semigroup (1.1) - (1-7) and a locally 
stable numerical method (2.1) - (2.2) consistent with it. Then the 
numerical method is convergent to the nonlinear semigroup. 

Proof. 

Let be given any regular Z ' C Z. T G (0, oo), compact K C X, and 
e > 0. Further, let us take any t G [0,T], At > 0, n G N, with 
uAt <T, emdu e K r)X[nX 

We note that 

(3.12) Cl^u - E{t)u = {ClfU - E{nAt)u) + {E{nAt)u - E{t)u) 
while 

Cl^u - E{nAt)u = 

(3.13) 

= j:i<p<n{Cl7CME{{p - l)At)u - Cl;'E{At)E{ip - l)At)u} 
However, the set 

(3.14) K" ^{CasE(s)v I s,As e[0,T], V e KnX'J 

is compact, since it is the image through the continuous mapping (2.1) 

- (2.2) of [0, T] X K' which in view of (3.7) is compact. Furthermore, 
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in view of (2.2), we have 
(3.15) K'CK" 

And now, the assumption of local stability applied to Z', T and K 
yields L',p' > 0, such that for each 1 < p < n, we have 

\\Cl-''CAtE{{p-l)At)u-Cl-PE{At)E{{p-l)At)u\\ < 

(3.16) 

< L'\\CAtE{{p - l)At)u - E{At)E{ip - l)At)u\ \ 

provided that 

(3.17) CME{{p-l)At)u, E{At)u, E{{p 
and 

(3.18) \\CAtE{{p - l)At)u - E{At)E{{p 

Here, the fact that (3.17) is satisfied follows easily from the assump- 
tion that u e X nAtJ ^^^^ (3.15) and (2.4), (2.5). 
As for (3.18), we note the following. The assumption of consistency in 
(3.11) applied to Z', T, K and e yields 5 > 0, such that for 1 < p < n, 
we have 

(3.19) \\CAtEiip - l)At)u - E{At)E{{p - l)At)u\ \ < eAt 
as soon ets At < 5. However, we can obviously assume that 

(3.20) At<d =^ eAt<p' 
and then (3.18) holds. 

Now (3.13), ( 3.16) and (3.19) give 

(3.21) At<5 ^ \\CltU - E{nAt)u\ \ < nL 'eAt < L 'Te 

while the assumption of well posedness applied to the compact Z ' fl 
([0, T] X K) results in the uniform continuity of the corresponding re- 



-i)Ai)Kei^"nx;„_^)^, 
- l)At)u\ \ < p' 
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striction of the mapping (1.6), thus leading to 5' > 0, such that 
(3.22) \t-nAt\<5' ^ \\E{nAt)u - E{u)\\ < e 
Finally, (3.21), (3.22) give 
WC^tU - E{t)u\\ < {l + LT)e 



as soon as 

At<5, \t-nAt\<5' 
and the proof of (3.1) is completed. 

Remcirk 3.1. 



□ 



As shown in Rosinger [2, pp. 67-69], Theorems 3.1. and 3.2. above 
contain as a particular case the Lax-Richtmyer Linear Theorem in its 
somewhat stronger form as far as the requirements involved in that 
result are concerned. 



4. Conditions for Nonlinear Equivalence between 
Convergence and Stability 

The above Theorems 3.1. and 3.2. show that, under the respective 
conditions, we have the implications 

"local stabihty convergence distant stabihty" 

therefore, in view of Proposition 2.1., it follows that 

"stability =^ convergence" 

In this way, in order to obtain an equivalence result between conver- 
gence and stability, that is, a stability type characterization of con- 
vergence, it suffices to find conditions - denoted by (C) in the sequel 
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- such that for numerical methods (2.1) - (2.2), the following implica- 
tion holds 

(*) CAt satisfies (C) \ 
=^ ( Ca* locally stable ) 
* *) C'a* distantly stable J 

And then, one easily obtains 

Theorem 4.1. ( Nonlinectr Equivalence) 

Given a well posed nonlinear semigroup (1.1) - (1-7) and a numerical 
method (2.1) - (2.2) which is consistent with it and also satisfies a 
condition of type (C). Then the numerical method is convergent to 
the nonlinear semigroup, if and only if it is stable. 

□ 

Remark 4.1. 

The conditions of type (C) in (4.1) are obviously bridging the gap 
between distant stability and local stability. And as seen above, such 
a gap occurs regardless of what seems as a rather natural approach 
to defining the concept of stability of a numerical method associated 
with a nonlinear semigroup. 

□ 

Here we shall present an explicit formulation of one specific condition 
of type (C) . Needless to say, there may be many other such formula- 
tions, some of them more useful in applications than other ones. 

In condition (2.8) defining distant stability, let us, for given Z',T,K, 
and p denote by 

(4.2) L"{Z',T,K,p) 

the smallest L " for which that relation holds. Then obviously L "{Z T, K, p) 
is a decreasing function of p, for Z' ,T and K fixed. Therefore, one 
obtains easily 
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Proposition 4.1. 

Given a numerical method (2.1) - (2.2) which is distantly stable. Then 
this numerical method is also locally stable, if and only if 

V regular Z' QZ, T e (0, oo), compact K QX : 

(4.3) 

\m^p^^L"{Z',T,K,p) < oo 

□ 

In view of (4.1), we obtain 

Corollary 4.1. (A condition of type (C) ) 

The above condition (4.3) is a condition of type (C) for numerical 
methods (2.1) - (2.2). 

5. Continuity and Compactness, instead of Completeness 
and Boundedness Are Relevant to Convergence 
and Stability 

Given a well posed nonhnear semigroup E{t) in (1.1) - (1.7), and a 
numerical method Ca* in (1.2), (2.2) consistent with it. 

Let us start by recalling in the particular linear case, the case in which 
the Lax-Richtmyer Equivalence Theorem alone operates, the respec- 
tive linear stability condition on the linear numerical method, valid on 
a given time interval [0, T] , namely 

3 L>0 : 

(5.1) V At>0, neN : 

nAt<T =^ \\Clt\\<L 
This simply means that the infinite family of linear operators 

(5.2) ( CXt I At > 0, n e N, nAt <T ) 
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has to be uniformly bounded on the normed space {X, \ \ ||). 

Second, and as mentioned above, see also Rosinger [1-9], the imphca- 
tion "convergent =^ stable" in the Lax-Richtmyer Linear Equivalence 
Theorem makes essential use of the Principle of Uniform Boundeed- 

ness of Linear Operators on Banach Spaces. Thus, as is well known, 
that theorem cannot operate in arbitrary normed spaces {X, \ \ ||), but 
only on Banach spaces, that is, normed spaces which are complete. 

This is, therefore, the way boundedness and completeness are essen- 
tially involved in the Lax-Richtmyer Linear Equivalence Theorem. 
And as mentioned, see also Rosinger [1-9], that fact leads to major 
troubles when nonlinear extensions of that theorem are attempted. 

On the other hand, a simple look at the phenomena involved in the 
approximation of semigroups by numerical methods does immediately 
reveal that, instead of boundedness and completeness, two other topo- 
logical features are in fact essentially involved, namely, continuity and 
compactness. And that fact, as seen above, as well as in Rosinger [1- 
9] , can help in finding nonlinear generalizations of the Lax-Richtmyer 
Linear Equivalence Theorem. 

Let us now present some of the main related details. 

As is well known, a vast and particularly important case of approxi- 
mation of semigroups by numerical methods occurs when the respec- 
tive semigroups correspond to classical solutions of nonlinear evolution 
PDEs. And typically in such cases, the respective PDEs arc known to 
have exact solutions, since otherwise, their numerical approximation 
would of course be pointless. Furthermore, a most important case of 
such known to exist exact solutions is when, in addition, they are also 
known to be classical, that is, sufficiently smooth. 

Consequently, the spaces X in which the respective nonlinear semi- 
groups operate, see (1.1) 

(5.3) E{t) -.Xt — > X, te [0, oo) 
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can be assumed as given by spaces of sufficiently smooth functions 
on suitable Euclidean domains, the domains on which the respective 
nonlinear evolution PDEs are defined. Thus, when endowed with typ- 
ical norms, such spaces X turn out not to be complete. Furthermore, 
by completing these spaces X in such norms, one enlarges them con- 
siderably, by adding to them large sets of highly nonsmooth functions. 

On the other hand, having assumed the existence of classical, that 
is, smooth enough solutions which are elements of such spaces X of 
smooth enough functions, there is simply no need to further extend 
those spaces, and in particular, there is no need to complete them in 
any way. 

And in fact, even in the general case when the nonlinear semigroups 
would not necessarily be associated with evolution PDEs, their defi- 
nition in (1.1) - (1-7) does not in any way require the completeness of 
the normed spaces {X, \ \ ||) on which they are defined. 

There has also been a rather different line of argument which is claimed 
to motivate the need for dealing with complete, that is, Banach spaces 
{X, II II). Namely, it is claimed that the inevitable presence of round- 
off errors when implementing numerical methods on digital computers 
leads to that necessity. 

As it happens, however, such arguments prove to be fallacious, see 
Rosinger [2, pp. 231-242], [4,5,7,9]. Indeed, the way which is typically 
suggested for the inclusion of the essentially nonlinear effects of prop- 
agation of round-off errors is so unrealistic that it leads to strictly 
better results in the case of numerical solutions of PDEs, than the 
well known best possible ones in the case of the numerical solutions 
of ODEs, see Rosinger [2,4,5,7,9], Isaacson- Keller. 

So much, therefore, for the requirement of competeness in the general 
nonlinear situation, that is, for the normed spaces {X, \ \ ||) to be in 
fact Banach spaces, a requirement which gets imposed only by the 
way the implication "convergent =^ stable" is proved in the particu- 
lar linear case of the Lax-Richtmyer Equivalence Theorem. 
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As for continuity, this condition in one or another form is typically re- 
quired both on the nonlinear semigroups and the numerical methods 
aimed to approximate them, as seen in (1.6) and (2.1). 

And once such a continuity is given, the presence of compactness is im- 
mediate. Indeed, for every given time interval [0, T] and initial value 
u e Xt, the respective trajectory 

(5.4) {E{t) I te[0,T]} 

is always a compact subset of X, being the image of the compact in- 
terval [0,T] through the continuous mapping E{.)u. 

Therefore, the approximation of the nonlinear semigroup E{t) by the 
numerical method Ca* is in fact but the approximation of the compact 
set (5.4). 

By eliminating completeness, one avoids having to deal with large ad- 
ditional sets of elements in the normed spaces (X, || ||), this being a 
considerable advantage. 

Second, in infinite dimensional normed spaces, typical when the non- 
linear semigroups are associated with evolution nonlinear PDEs, com- 
pact sets are bounded, but not the other way round. Furthermore, 
compact sets are far smaller than bounded sets in general. Conse- 
quently, formulating stability conditions involving compact, and not 
bounded sets, lead to far weaker such conditions. 

In conclusion, there is a whole variety of advantages which result from 
formulating equivalence theories between stability and convergence in 
terms of continuity and compactness, and not in terms of complete- 
ness and boundedness, see Rosinger [1-9] . 
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